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Abstract 

If the gravitino mass is in the region from a few GeV to a few lO's GeV, the scalar lepton X such 
as stau is most hkely the next hghtest supersymmetry particle. The negatively charged and long- 
lived X~ may form a Coulomb bound state {AX^ ) with a nucleus A and may affect the big-bang 
nucleosynthesis through catalyzed fusion process. We calculate a production cross section of ^Li 
from the catalyzed fusion (^HeX~) -|- d — > ^Li -|- X~ by solving the Schrodinger equation exactly 
for three-body system of ^He, d and X. We utilize the state-of-the-art coupled-channel method, 
which is known to be very accurate to describe other three-body systems in nuclear and atomic 
reactions. The importance of the use of appropriate nuclear potential and the exact treatment 
of the quantum tunneling in the fusion process are emphasized. We find that the astrophysical 
5-factor at the Gamow peak corresponding to T = 10 keV is 0.038 MeV barn. This leads to the ^Li 
abundance from the catalyzed process as ^Li|cBBN — 4.3 x 10~-'^^(D/2.8 x 10^^)([ri^_ /s]/10^^^) in 
the limit of long lifetime of X. Particle physics implication of this result is also discussed. 
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I. INTRODUCTION 

The gravitino is the most important prediction of supergravity (SUGRA) Its mass 
is expected to be in a wide range, 1 eV— 100 TeV, depending on the mediation mechanism 
of supersymmetry (SUSY) breaking effects. In particular, the gravitino with a mass of 
0(1) MeV— 0(10) GeV is very interesting in the sense that it is most likely the stable and 
lightest SUSY particle (LSP) and can be a candidate of cold dark matter in the universe. 
Moreover, if the gravitino mass lies between a few GeV and a few lO's GeV and the next 
LSP (NLSP) is a charged scalar lepton (X) such as stau (f), the SUGRA can be tested in 
high-energy collider experiments 2j. The key observation is that the gravitino mass ^3/2, 
the stau mass and the stau lifetime could be measured in high precision if m3/2 
is of 0(10) GeV (cf. 31). Then the Planck scale MpL is extracted through the relation, 
^^xp.) _ x^(77j(^^p )^ Mpl). This independent determination of MpL will provide us 

with a crucial test of the SUGRA. 

For the gravitino mass in such an interesting region as discussed above, the stau has 
necessarily a long lifetime. For example, it is of order one month for m3/2 — 10 GeV and 
nix — 100 GeV. However, such a long-lived particle is potentially dangerous in cosmology. 
If it decays after the big-bang nucleosynthesis (BBN), the decay products destroy the light 
elements and ruin the success of the SBBN (standard BBN) {4, 5|. More seriously, if the long- 
lived particle is negatively charged as X~ , it forms bound states together with positively 
charged nuclei, which leads to an enhancement of some nuclear reaction rates. That is, 
X~ plays as a catalyzer of nuclear fusion (catalyzed BBN or CBBN in short) and affects 
cosmological nuclear abundances.^ 

In a recent article Pospelov has argued that too much ^Li is produced through the 
CBBN if the lifetime of X~ is long enough, > 10^ — 10*^ sec, and its abundance is large 
enough, n^_/s > 10^^'', where n^_ and s are the number density of X~ and the entropy 
density of the universe, respectively.^ The production rate of ^Li estimated in Ref. [3] is 
based on a naive comparison of the standard process ^He+c/ ^Li+7 and the new catalyzed 
process (^HeX~) +d — ^ ^Li + X^ where (^HeX~) is the Is Coulomb bound state of ^He and 



^ Massive and long-lived X produced in accelerators may be of practical use as a catalyzer for D-D and 

D-T fusions similar to the muon catalyzed fusion [gj. 
^ Other effects of the X~ bound states have been considered in 
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X~ . It is assumed that the standard process is dominated by the E2 transition^ induced 
by the interaction QijViEj with Qij being the quadrapole operator. Then, by applying the 
same interaction to the photonless process with Ej replaced by the Coulomb field associated 
with X~ , it was concluded that the astrophysical S-factor of the new process at zero incident 
energy is about 0.3 MeV barn. However, the assumptions adopted in the above estimate do 
not have firm ground: For example, the angular momentum of the initial ^He-rf system in the 
catalyzed process is dominated by zero, L = 0, while the angular momentum of the initial 
system in the standard process is dominated by L = 2 (L = 1) for the E2 (El) transition. 
Such a kinematical difference invalidates the use of the SBBN process to estimate the CBBN 
process. Furthermore, the quantum interplay between the ^He — d nuclear fusion and the 
tunneling of d through the Coulomb barrier plays an important role in low-energy catalyzed 
fusion and cannot be treated in a perturbative manner. 

The purpose of this paper is to solve the Shrodinger equation for the three-body system 
(^He, d and X~) exactly and derive a reliable S'-factor for the CBBN. The method 



we adopt is the state-of-the-art coupled-channel technique deve 
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oped by two of the present 



14| . Since we can treat the 



authors (M.K. and Y.K.) together with E. Hiyama jll . 
catalyzed process directly, we do not need to refer to any of the SBBN processes. Also, the 
method has been already proven to be highly accurate and useful in atomic and nuclear 



physics (reviewed in ij]). We find that the obtained astrophysical S'-factor at E = 36.4 



keV (the position of the Gamow peak for T = 10 keV which is relevant to the CBBN) 
is 0.038 MeV barn. This is about 10 times smaller than the estimate in Ref. j?] at the 
same energy. For long lived stau, our S-factor leads to the ^Li abundance from CBBN as 
•^LilcBBN - 4.3 X 10"^^(D/2.8 x 10'^){[n^_ / s]/10^^^). Therefore, an observational upper 
bound on the ®Li abundance ^Li < 6. 1 x 10"^^ (2cr) sl] leads to a bound on the X~ abundance, 
^x- / s < 1-4 X 10^^^. This requires a dilution of the relic stau by some entropy production 
at late time by a factor of A ^ (300 — 600) x (m~/100 GeV). Such a dilution factor can 
easily be consistent with the thermal leptogenesis |l5| for a reheating temperature Tr > 10^^ 
GeV. 

The organization of this paper is as follows. In Sec. II, we summarize the basic reaction of 



^ This particular assumption itself may not be justified since Ef transition could be comparable or even 
larger than E2 at low energies. See e.g. [lo| . 
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the X~ catalyzed fusion as well as its atomic analogue, the muon transfer reaction. In Sec. Ill, 
we present our method of solving the ^He-ci-X" three-body problem. In Sec. IV, we show 
our result of astrophysical ^-factor as a function of the incident energy. Sec.V is devoted 
to conclusion and discussions in which the ^Li abundance in CBBN and its implication to 
particle physics are also mentioned. 

II. THE X- CATALYZED PROCESS 

As we have mentioned in the Introduction, we will treat the following process directly in 
a fully quantum mechanical manner: 

(^HeX-) + rf^ ^Li + X- + l.lMeV. (2.1) 

As shown later in Section V, the relevant temperature T for CBBN is about 10 keV. This 
corresponds to the Gamow peak of the incident channel located at the cm. energy =36.4 
keV with the full 1/e width of 44 keV. Thus the cm. energy of the incident channel will 
be typically up to 100 keV. We adopt the ^He — d cluster model in which the ^Li nucleus 
is a bound state of a ^He nucleus and a deuteron. This model has extensively been utilized 
in the studies of the structure and reactions of light nuclei and is very well established (see 



e.g 
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17(1). Such a cluster-model treatment allows us to investigate the reaction (12. ip as 
a quantum three-body problem, ^He + d + . The rearrangement of the three ingredients 
takes place during the reaction process through the long range Coulomb interaction and the 
short range nuclear interaction. 



Although solving the three-body problem accurate 
authors (M.K. and Y.K.) have an experience [l2, [l3 



y is quite an elaborate task, two of the 
Q to have solved a similar subject. 



the muon transfer reaction, 

(rf/i")i, + t ^ (t/i")i3 + rf + 48 eV (2.2) 

at incident energies of 0.001 — 100 eV in the context of the muon catalyzed D-T fusion cycle 
jisl . Among the calculations of this reaction in the literatures, their three-body coupled- 
channel method was found to provide with the most accurate result (see §8 of [14']). 

In the present paper, the same coupled-channel method is applied to the reaction, Eq. 
(12. ip . We note here that it is very important to employ an appropriate nuclear interaction 
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between "^He and deuteron in treating the reaction; we choose the interaction which can 
reproduce the binding energy and the charge form-factor of ^Li as well as the low-energy 
"^He — d scattering phase shift. 

III. COUPLED CHANNEL METHOD 

Following Ref ■ . we briefly explain the three-body coupled-channel method to investi- 
gate the reaction (12.11) . As shown in FigUl we consider all the Jacobi coordinates (rc, Rc) 
of the three possible sets, c = 1,2 and 3. The entrance channel (^HeX~) + d and the exit 
channel ^Li + X~ in our reaction are best described by the coordinates for c = 1 and those 
for c = 2, respectively. Since we are interested in the incident energy of the entrance channel 
below 100 keV, there is no other open channel than the above two. All the excited states 
of (^HeX~) and ^Li as well as all the states of {dX~) can be excited only virtually in the 
intermediate stage of the reaction. To describe such an intermediate state where particle- 
rearrangement takes place, it is convenient to utilize the coordinates for c = 3 together with 
those for c = 1, 2. 



4 4 4 

He He He 




c=1 c=2 c=3 

FIG. 1: Three sets of Jacobi coordinates of the ^He + d + X^ system. The entrance (exit) channel 
is described by the coordinate system of c = 1 (c = 2). 
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A. Three-body Schrodinger equation 



The Schrodinger equation for the total wave function \I/jm of the ^He + d + X system 
having an angular momentum J and its z-component M is given by 

(i/-^tot)^JM = 0, (3.1) 

with the Hamiltonian, 

^ = " 2^^^^= + V^^He-x(ri) + V.nUr2) + V,.x{r,). (3.2) 

As far as we use the reduced masses (m^ and associated with the coordinates (r^ and 
Rc), every choice of c in the kinetic term is equivalent. Va^-bX''^c) denotes the potential 
between the particles Ac and Be {Ai-Bi = ^He-X, A2-B2 = '^He-rf, A3-B3 = d-X) and will 
explicitly be given below. 

Spin of the deuteron is neglected, and therefore the angular momentum of the ground 
states of ®Li is zero as well as that of (^HeX~). We denote the ground-state wave function 
of (^HeX~) in c = 1 by (ri) and its eigenenergy by e^^l , and similarly, that of ®Li in 
c = 2 by 0g^s.(r2) and its eigenenergy by e^^} . They are obtained by solving 



2mr 



Vrl + Va^-bM - 41 ] = 0. (c = 1, 2) (3.3) 



Center-of-mass scattering energy of the channel c associated with the coordinate Rc, say 
i?c, is introduced by Ec = -Etot — ^^gl together with the corresponding wave number kc by 
h^k1/2Mc = Ec{c = 1,2). For a given total energy i^tot, the Schrodinger equation (13. ip 
should be solved under the scattering boundary condition: 

hm = (pii{rM'\kMSci- Ms^^cU^/\kcR^^^ (c=l,2) (3.4) 

where u^j^\kR){= {Gj{kR) ± iFj{kR))/kR) are the asymptotic outgoing and incident 
Coulomb wave functions. Sf^^ is the S-matrix for the transition from the channel 1 to 
c and Vc is the velocity of the channel c. By introducing a simplified notation, E = Ei and 
k = ki, the cross section of the rearrangement process (12. ip is given by 

00 

^i^2{E) = -2T.i'^J+l)\SL2\\ (3.5) 

J=0 
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and the astrophysical S'-factor is derived from 

SiE) = a,^2{E)Eexp{27rriiE)), (3.6) 

where r]{E) is the Sommerfeld parameter of the entrance channel. 

The three-body wave function which describes the transfer reaction (12.11) and the elastic 
(^HeX^) + d scattering simultaneously is written as 

^JM = 02.(ri) x%{Ri) + 0g.(r2) xS/(R2) + ^itr^ . (3.7) 

The first and the second terms represent the open channels, c = 1 for (^HeX^) + d and 
c = 2 for ^Li + X~. The factors Xjm(R'c)(= Xj^(-Rc) ^jj\/(R-c) ) describe the scattering 
waves along the coordinates Rc and are to be solved under the boundary condition (13.41) . 
The third term, stands for all the closed (virtually-excited) channels in the energy 

range of this work; in other words, this term is introduced to represent all the asymptotically- 
vanishing three-body amplitudes that are not included in the first two scattering terms. For 
example, the third term describes such an effect that the incoming deuteron attracts the 
^He in the Is-orbit around X~ and distorts the orbit before picking up the "^He. 

Since vanishes asymptotically, it is reasonable and useful 12,[l3|, [l4| to expand it 

in terms of a complete set of L^-integrable three-body basis functions, {^jm,u', v = 1 — t'max}, 
spanned in a finite spatial region (see Sec.3.C): 

^fM"'^=Ebj.^JM,.. (3.8) 

Equations for x^j\Ri), xf \R2) and the coefficient bji^ are given by the z/jnax + 2 simultaneous 
equations 

(</.(l(rc)>JM(Rc) |i/-i?tot|^JM)r,,R^ = 0, (c=l,2) (3.9) 

and 

( $ JM, u\H-E^,,\ ^JM ) = 0. (Z/ = 1 - Z/,„ax) (3. 10) 

Here, ( )^ g denotes the integration over and Rc. 

Since $jm, u are constructed so as to diagonalize the three-body Hamiltonian as 

($jAf,z.| H \^JM,u') = Ejy5yu', (Z^, y' = 1- Z/max) (3.11) 

the coefficients hji, can be written, from Eqs. fl3.10p . as 

hj, = ^ {^JM,. I H - i?tot I + 0S.X!/!). (^ = 1 - ^max) (3.12) 



Eju — E, 



tot 
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Inserting Eqs. fl3.12p into bji, in \I/jm in Eqs. fl3.9p . we reach two coupled integro-differential 
equations for Xj\Ri) and x^j\R2)- 



The integro-differential equations, though not recapitulated here (c.f. §8 of [14| for them), 
are solved by using both the direct numerical method (finite-difference method) and the 
Kohn-type variational method |14J . |19I | , and we have obtained the same result for the incident 



energies relevant to CBBN, > 10 keV. The coupling between the entrance and exit 
channels as well as the contribution from the closed channels ^'^f^^^'^^ were found to be 
significantly large as will be discussed later. 



B. Nuclear potentials 

It is essential for the three-body calculation to employ appropriate nuclear interaction 
between ^He and deuteron which governs the ^He-transfer process. In this subsection, we 
define the potentials V4He-x('"i), ^He-d('^2) and Vd-xij^) in our Hamiltonian (13. 2p . First, we 
assume Gaussian-shape charge distributions of ^He and deuteron as 2e{T(h\)~^^'^e~^'^^^^^'^ and 
e{'nh1)^^/'^e~''^/^^^'^ ^ respectively. We take hi = 1.37 fm and 63 = 1.75 fm, which reproduce 
observed r.m.s. charge radii, 1.68 fm of "^He and 2.14 fm of deuteron. The potential between 
^He and X~ is then given by 

l/.H._v(n) = -2e=5^1<^, (3.13) 

ri 

and that between deuteron and X~ is written as 

V,,(n,) = -e^^^i<r5M. (3.14) 

where erf(a;) = ^ Jq e~*^dt is the error function. Energy of the (^HeX~)is state is e^^l = 
—337.33 keV and the r.m.s. radius is {rf )^/^ = 6.84 fm. 

The potential V4He-d('"2) is a sum of the nuclear potential, Kne-dl^s), and the Coulomb 
potential, KHe-d(^2)- The latter is given by 



, eTf(r2/Jbj + bl) 
V.S,Ur2) = 2e' '\ (3.15) 

The nuclear potential is assumed to have a two-range Gaussian shape as 

V.l^Ar2) = vo e-^''/''^" + v', e'^''^/'^')' (3.16) 
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with a = 0.9 fm, vq = 500.0 MeV, a' = 2.0 fm and v'^ = -64.06 MeV. The first term, 



a repulsive core, is introduced to simulate the Pauli exclusion principle that nucleons in 
the incoming deuteron should not occupy the nucleon s-orbit in the ^He nucleus during the 



reaction process (see e.g. 



16| for this role of Pauli principle) 





e'^' [MeV] 



FIG. 2: Left panel: Charge form factor of the electron scattering from ^Li. The calculated values 
(experimental data 20|]) are shown by the solid line (filled circles). Right panel: The s-wave phase 



shift 6q of the ^He + d scattering at cm. energy e^m. < 10 MeV. The calculated values are shown 
by the solid line, while the data from the phase-shift analysis are shown by filled circles [21] and 



by open circles 



The above parameters of V^^^_^{r2) were so determined that the solution to the 
Schrodinger equation (13.31) could reproduce observed values of the energy e^}^ = —1.474 
MeV and the r.m.s. charge radius 2.54 fm of the ground state of ^Li. Furthermore, the 
charge density of ^Li reproduces observed charge form factor of the electron scattering from 
^Li as shown in the left panel of Figj2l Simultaneously, use of the potential V4He-(i(?"2) ex- 
plains the low-energy s-wave phase shifts of the ^He+c? scattering as shown in the right panel 
of Figj2l Here, it should be noted that simple attractive potential without a repulsive core 
leads to the phase shift which increases with increasing E and cannot explain the observed 
data. We thus have a good '^He — d potential V4He-d(?"2) in order to perform a precise study 
of the three-body reaction (12.11) . The potential V4He-(i('"2) and the wave function 0^^] (r2) are 
illustrated in Figl3] in which r2 is denoted by r for simplicity. 
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FIG. 3: The potential V4jjg_^(r) between "^He and deuteron (the soUd Une). To show its Coulomb 
barrier tail, the same potential scaled by 100 is shown by the dotted line together with the typical 
incident kinetic energy Eq x 100 denoted by the arrow. The dashed line shows the ^He — d relative 
wave function in the ^Li ground state in an arbitrary unit, 0g^l(r). 

C. Three-body basis functions 

The L^-integrable three-body basis functions {^jMy, u = \ — z/max} used in f l3.8l) to 
expand ^^f^^'^'^^ are introduced as follows [l^: $jm, are written as a sum of the component 
functions in the Jacobi-coordinate sets c = 1 — 3 (FigJT]), 

^jM,u = ^%Jv^,R^) + $S,,.(r2,R2) + $S,,,(r3,R3) . (3.17) 

Each component is expanded in terms of the Gaussian basis functions of the coordinates Tc 
and He'. 

"^JM.i/l^c, R-c) = X! '^Ju,nJc,NcLc 
nJcNcLc 

where 

0^,„(r) = r' e-(^/'^")' r,™(?), r„ = n a"-\ {n = 1 - n^ax) (3.19) 

^W(R) = e-(^/^^)' Ylm(R), Rn = Ri (iV = 1 - iV^,J (3.20) 

The Gaussian ranges are postulated to lie in geometric progression. Basis functions so 
chosen are suited for describing both short range correlations (mainly due to the ^He — d 
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G 



JM 



(c=l-3), (3.1^ 



nuclear interaction) and long range asymptotic behavior simultaneously [IJ], and therefore 
they are efficient to describe the three-body configurations (closed- channel contribution) in 



the interaction region in the intermediate stage of reactions 
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Q, Q. The coefficients 



^jl njc NcLc fl3.18p are determined by diagonalizing the three-body Hamiltonian H as 
Eqs.dSH]). 

In the calculation for J = 0, we took Ic = Lc = 0, 1, 2 and nmax = A^max = 15 for c = 1 — 3. 
Therefore, total number of the three-body Gaussian basis functions (f)'^j^{rc) tp'^^i^(Rc) 
to construct {^jm,u} amounts to u^ax = 2025, which was found to be large enough for the 
present calculation. As for the Gaussian ranges, we took ri,rn^^^, Ri and -RA^max to be 0.5, 
15.0, 1.0, 40.0 fm, which are sufficiently precise for the present purpose. The expansion (13. 8p 
was found to converge quickly with increasing u, and u < 100 {Ej^ < 1 MeV above the 
(^HeX~) — d threshold) is very sufficient. 



IV. RESULT OF THE REACTION RATE 

Table [I lists the cross section cxi^2{E) and the astrophysical S'-factor S{E) of the reaction 
(12. ip obtained by the full coupled-channel calculation at several energies around Eq = 36.4 
keV. In FiglU S{E) is shown together with the Gamow peak. The S'-factor increases with 
decreasing energy E since there are three-body bound states below the ^Li-X~ threshold. 

n 

Note that S{Eg) estimated in [7] is an order of magnitude larger than the value given by 
the present exact calculation. 

Since the S'-factor can be approximated around Eq as 

S{E)c:^S{EG) + a{E-EG), (4.1) 

with a = {dS /dE)^^, the reaction rate, using Eqs.(4.56) and (4.74) of [2^, is expressed as 

NA{av) = 6.24x10^ S (Eg) {1-0.0718 aTg/ S{Eg))Tq'^^'^ exp {-5.33Tq^^'^) cm^s~^mor\ 

(4.2) 

where Na is the Avogadro constant and S{Eq), a and Tg are to be given in units of MeVbarn, 
barn and 10^ K, respectively. Taking S{Eq) = 0.0380 MeV barn and a = —0.18 barn, we 
obtain 

Na{(Jv) =2.37x 10^(1- 0.34 T9)T9"^/^exp (-5.33X9"^/^) cm^s'^moPV (4.3) 
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TABLE I: The cross section ai^2{E) and the astrophysical S'-factor S{E) of the reaction (^HeX ) + 
d — 5> ^Li + X~ obtained by the fuh coupled-channel calculation. The most effective energy Eq for 
T =10 keV is 36.4 keV. 

E [keV] cri^2 [barn] S [MeV barn] 



10 


3.85 X 10' 


-6 


0.0426 


20 


1.09 X 10" 


-4 


0.0410 


36.4 


6.88 X 10- 


-4 


0.0380 


50 


1.41 X 10" 


-3 


0.0357 


100 


3.50 X 10- 


-3 


0.0286 




50 100 

E [keV] 

FIG. 4: The astrophysical 5-factor (the solid line) obtained by the full coupled-channel calculation. 
The dashed line illustrates the Gamow peak (in arbitrary units) for T = 10 keV with the peak 
maximum at Eq = 36.4 keV. 

To check the necessity of the exact coupled- channel treatment (13.71) — fl3.12p . we have cal- 
culated the cross sections in the following approximations: (i) A coupled-channel treatment 
with only the entrance and exit channels where the closed channel '^^jm"^'^^ in Eq. ( 13.70 is ne- 
glected, and (ii) an one-step approximation where the transition from the entrance channel 
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to the exit channel in the case (i) is treated in the leading order perturbation. For the case 
(i), the cross section o"i^2 become approximately 3 times smaller than those in Table [H For 
the case (ii), the cross section becomes an order of magnitude larger than that in the case 
(i). Therefore, these approximations are not justified. 

We have also checked that the cross section cri_»2 is dominated by the total orbital angular 
momentum J = 0, and contribution from the higher angular momenta J > 1 is three orders 
of magnitude smaller. Furthermore, any pair among ^He, d and X~ has dominantly zero 
orbital angular momentum, i.e. s-wave. This is in sharp contrast to the SBBN reaction 
'^He + d ^ ^Li + 7 where the entrance channel is dominated either by d-wave (the E2 
transition) or by p-wave (the El transition) from the selection rule. 



V. CONCLUSIONS AND DISCUSSION 



Since we consider the gravitino masses of order 10 GeV, the scalar lepton X (such as 
stau) is most likely the next lightest supersymmetry particle and could have a long lifetime 
more than 10^ sec. Thus, the negatively charged X~ forms a bound states {AX~) with 
positively charged nuclei A during the big-bang nucleosynthesis (BBN). Then, some nuclear 
abundances are enhanced through the catalyzed process. 

In this Letter, we have calculated a production cross section of ^Li from the catalyzed 
process (^HeX~) + d ^ ^Li + X~ by exactly solving the Schrodinger equation for three- 
body system of ^He, d and X. We utilize the state-of-the-art coupled-channel method, 
which is known to be very accurate to describe other three-body systems in nuclear and 
atomic reactions. The use of appropriate nuclear potential and the exact treatment of the 
quantum tunneling in the fusion process turned out to be important. We have found that 
the astrophysical S'-factor at the Gamow peak corresponding to T = 10 keV is 0.038 MeV 
barn as shown in Table [U Fig. [Hand Eq. (14. 3p . They are the main results of this Letter. 

Before closing, let us briefly discuss the particle physics implication of our result. With 
the cross section (erf) in Eq. (14.31) . the ^Li production via CBBN is described by 

d 



dt 



= nBs{crv)'D, (5.1) 

CBBN 



where D = Ud/nB and ^Li = n^uJiiB are the abundances of these elements normalized by 
the baryon number density. The ratio of the number density of the bound states tt-bs to the 
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FIG. 5: (a) The number fraction of particles which form bound states with ^He as a function 
of temperature T. tx = (1,2,4,00) x 10'^ sec from bottom to top. (b) The contour plot of 
^Li abundance produced by the CBBN process in {tx, ^x-fs) plane. The solid lines represent 
^Li = 10^^^, 10^^^, and 10^^^ from the bottom to the top. The dashed lines represent the thermal 
relic abundance of stau for m~ = 100 GeV (below) and 300 GeV (above). 



entropy density is given by ^ 

^Bs nx- 



eM-t/rx) (52) 



s s initial 1 + n-i(mQ,T/27r)3/2exp(-i?fe/T)' 

where tIq, and rria are the number density and the mass of the ^He, and Eb GeV is the 
binding energy of (^HeX^). The fraction of X~ which form bound states with ^He is shown 
in Fig. [5](a), as a function of temperature T for various lifetimes of X.^ The bound state 
abundance is peaked at around T = 10 keV for > 1000 sec,^ which justifies the expansion 
in Eq.dlT]). Thus, the CBBN occurs at t > 1000 sec, when the SBBN processes (t < 1000 
sec) are already frozen. At such late time, the effect of the dissociation ^Li + p He +^ He 

n 

is also negligible [Sj, and hence we have neglected it in Eq. (15. ip . 

Fig. [5t^b) shows the contour plot of ^Li abundance produced by the CBBN process, as a 
function of initial X~ abundance nx-/s and the lifetime t^. In the limit of long lifetime, 



^ We have assumed that the reaction "^Hc + X <-> (^HeX ) is in chemical equiUbrium. For more detailed 

analysis of the evolution of the bound state abundance, see Ref. Q. 
^ For our numerical calculation, we use Y — Aua/uB — 0.25, D = 2.78 x 10~^, ub/s = 0.87 x 10"^", 

s = 1.715r3, rria = 3.73 GeV 12^ and Eb = 337.33 GeV (see the previous section). 
^ For < 1000 sec the CBBN effect is negligible. See Fig.IDJb). 
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one obtains 

^Li 



_ 4.3xlO--f^V ^) . (5.3) 

CBBN I 10-16 / V2.78 X 10-5/ ^ ' 



Therefore, the observational upper bound on the ^Li abundance ^Li < 6.1 x 10^^^ (2(t) [sI 
leads to a bound on the X~ abundance, nx-/s < 1.4 x 10~^^. On the other hand, the 
thermal relic abundance of stau is given by n~- / s ~ (4 — 8) x 10-^''(m~/100 GeV) 
(Conservative cases are shown in Fig. [5]^b) as the dashed lines for m~ = 100 GeV and 300 
GeV). Therefore, in the limit of long life time, an entropy production with a dilution factor 
A ~ (300 — 600) X (m~/100 GeV) for the primordial stau abundance is necessary to avoid 
the overproduction of ^Li.^ 

On the other hand, this means that the primordial baryon asymmetry is also diluted 
by A. Standard thermal leptogenesis [l^, an attractive mechanism for baryogenesis, can 
normally explain the observed baryon asymmetry if the reheating temperature satisfies Tr > 



2 X 10^ GeV 27|. Since the maximal baryon asymmetry is proportional to Tr in thermal 
leptogenesis, this means that, for a late time dilution A ^ (300 — 600) x (m~/100 GeV), the 
reheating temperature of Tr > (0.6 — 1.2) x 10^^ GeV x (m~/100 GeV) is necessary. Such 
a reheating temperature can be obtained in inflationary models. 
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